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Abstract 

The Fredholm determinants of a special class of integral operators K supported 
on the union of m curve segments in the complex A-plane are shown to be the r- 
functions of an isomonodromic family of meromorphic covariant derivative operators 
T>\. These have regular singular points at the 2m endpoints of the curve segments 
and a singular point of Poincare index 1 at infinity. The rank r of the vector bundle 
over the Riemann sphere on which they act equals the number of distinct terms 
in the exponential sums entering in the numerator of the integral kernels. The 
deformation equations may be viewed as nonautonomous Hamiltonian systems on 
an auxiliary symplectic vector space M, whose Poisson quotient, under a parametric 
family of Hamiltonian group actions, is identified with a Poisson submanifold of the 
loop algebra gl R (r) with respect to the rational H-matrix structure. The matrix 
Riemann-Hilbert problem method is used to identify the auxiliary space M with the 
data defining the integral kernel of the resolvent operator at the endpoints of the 
curve segments. A second associated isomonodromic family of covariant derivative 
operators T> z is derived, having rank n=2m, and r finite regular singular points at the 
values of the exponents defining the kernel of K. This family is similarly embedded 
into the algebra gl R (n) through a dual parametric family of Poisson quotients of M.. 
The operators T> z are shown to be analogously associated to the integral operator K 
obtained from K through a Fourier-Laplace transform. 



1. Introduction 

la. Dual Isomonodromic Deformation Equations. 

The differential equations determining isomonodromic deformations of rational covariant 
derivative operators of the form 

V x = ^--Af(\), (1.1a) 
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3=1 A ^ 

5 = diag(A,.../? r ), Nj e [(r,C) 

were derived in [JMMS, JMU]. The residue matrices Nj(ai, (3 a ) may depend parametrically on 
the pole locations {a\, . . .a n } and the asymptotic eigenvalues {(3i, . . .(3 r }. This dependence 
must be such as to satisfy the commutativity conditions 

[V x ,V a .]=0, [V x , V Pa }=0, j = l,...n, a=l,...r, (1.2) 

where the differential operators V CCj , Vp a are defined by 

V aj := — + —^. 1.3a 
3 oaj A — a.j 

and E a is the elementary r x r matrix with elements 

(E a ) bc := 5 ab 5 ac . (1.4) 

In the Hamiltonian formulation [H] these relations, viewed as differential equations for the 
residue matrices Ni, are interpreted as a compatible system of nonautonomous Hamiltonian 
equations with respect to the Lie Poisson structure on the space (g[(r))* n = {N\, . . . N n } : em- 
bedded through formula (1.1b) into the loop algebra £)[r(?") with rational .R-matrix structure. 
They are generated by the Poisson commuting family of Hamiltonians {Hj, i^ }j=i...n,a=i...r 
defined by 

Hj^BN j )+±^^- (1.5a) 

i=i 3 i 



- (E; =1 ^) a6 (ELi^) 6 a 



-V \ ^ * 7 — -L J I , * 1% — J- ' OCL 

<*3 (^)aa + E — f^7J • ( L5b ) 

,• = 1 6 = 1 ^ a fJb 

The Poisson commutativity of the Hamiltonians follows from the .R-matrix structure on fll^(r) 
and implies that the differential 1-form on the parameter space {a t . . . a n , (3\, . . . (3 r } defined 
by 

n r 

u:=Y, H k da k + K adf3 a (1.6) 

k=l a=l 
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is closed 



du = 0. 



(1.7) 



This, in turn, implies the existence of a tau-function r(ai, . . . a n ; 0i, . . . f3 r ) in the sense of 
[JMU], determined, up to multiplication by a local constant, by 



uj = (i(lnr). 



(1.8) 



As shown in [H] it is natural to view the phase space as a Poisson quotient of a symplectic 
space M. = {F, G} consisting of canonically conjugate pairs (F, G) of rectangular N x r 
matrices, where 



N = ^ ki, hi := rankiVj 

i=i 

The quotient map is defined by the formula 



(1.9) 



Af(\) = B + F 1 (A-XI^G, 



(1.10) 



where A e gl(N,C) is a diagonal matrix with eigenvalues {oii\i=\... n and multiplicities {ki}. 
The Hamiltonians may be pulled back to M. through the projection map defined by (1.10), 
and hence determine corresponding nonautonomous Hamiltonian systems 



f) F f)G 

" {F, Hi}, ^L = { G ,H i }, i = l... 



dcti 
dF 

Wa 



{F, K a }, 



d(Xi 

dG 

Wa 



{G, K a } a = l,...r, 



(1.11a) 
(1.11b) 



whose integral curves project to solutions of the isomonodromic deformation equations (1.2). 

This leads quite naturally to a "dual" isomonodromic deformation system for the family 
of operators 



V Z :=§--M( Z ), 



M(z) := A + F(B - zI)- 1 G T = A + J2 "iV' 

a=l Z ^ 

M a e Q l(N,C). 



(1.12a) 
(1.12b) 



The corresponding deformations equations are 



V z , V c 



o, 



, j = l,...n, a = l,...r, 



(1.13) 
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where 

V ai -B-zEj-Y, E ' g °=' M ' ] E ' + E ' (£ °=' M ' ] E ' . (L14b) 

OCXi . =1 Oti Oij 

These have the same interpretation, as nonautonomous Hamiltonian equations on the phase 
space (gl(n))* r = {Mi, . . . M r }, generated by the Poisson commuting family of Hamiltonians 

K..^iAM.) + ±^^, (1.14a) 

t^f Pa- Pb 

Sj:=±0. (M a)j , + ± (£L ' ^ <ELl Mi) « . (1.14b) 

Pulled back to the space .M, these are shown in [H] to coincide with the Hamiltonians (1.5a), 
(1.5b), and hence the corresponding "dual" r-function r(«i, . . .a n ; (3\ . . .(3 r ), determined by 

n r 

d(lnr) =:= ^ H k da k + ]T K a d(3 a , (1.16) 

fc=l o=l 

coincides with r(cti, . . . a n ; fli . . . (3 r ). 

In the present work, using the approach of [IIKS] , it will be shown how a particular class 
of solutions to such isomonodromic deformation equations arises through the solution of an 
associated matrix Riemann-Hilbert problem (Theorems 2.1, 2.3), and how the corresponding 
r-functions may be identified with the Fredholm determinants of a simple class of "integrable" 
Fredholm integral operators (Theorem 2.5). We also show (Theorem 2.6) that the r-function 
obtained in this way may equivalently be interpreted as a r-function in the sense of Sato 
[Sa] and Segal and Wilson [SW]; that is, as the determinant of a suitably defined projection 
operator applied to the image of an element W of an infinite dimensional Grassmann manifold 
under the action of an abelian subgroup T r C <5l(r) of the loop group <25l(r). The "dual" system 
is then shown (Theorem 3.1), under suitable restrictions on the input data, to correspond to 
the integral operator obtained through a Fourier-Laplace transform. In this formulation, the 
element W of the Grassmannian is determined by the pole parameters {cti, . . . a n } of the first 
family of operators (1.1b), while the group element 7 G T r is parameterized by the asymptotic 
eigenvalues {f3±, . . . j3 r }. In the dual system, the roles of the element W and the group element 
7 of are interchanged. 
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lb. Integrable Fredholm Kernels and the Riemann—Hilbert Problem. 

Consider a p x p matrix Fredholm integral operator acting on C p -valued functions v(A), 

K(v)(A) = ^i^A,/z)v(/z)^, (1.17) 

defined along a piecewise smooth, oriented curve V in the complex plane (possibly extending 
to oo), with integral kernel of the special form 

A-(A,„) = q^M, (1.18) 

where f , g are rectangular r x p matrix valued functions, p < r. For the moment, we assume 
only that f and g are smooth functions along the connected components of V and, in order 
that K be nonsingular, we also require that the numerator vanish at least linearly in A — fi in 
the limit A — > fx, 

f T (A)g(A) = 0, (1.19) 
so that the diagonal values are given by taking the limit 

K(X, A) = f /T (A)g(A) = -f T (A)g'(A). (1.20) 

Fredholm determinants of operators of this type appear as generating functions for corre- 
lators in many integrable quantum field theory models [WMTB, JMMS, IIKS, IIKV, KBI] and 
as spectral distributions for random matrix ensembles [M, TW1-2, HTW], the most common 
case being p = 1, r = 2, with K a scalar integral operator. A first important observation to 
note [IIKS] is that the resolvent operator 

R:=(I-K) _1 K (1.21) 

is also in the same class. Thus, R may also be expressed as 

R(v)(A) = J R(X,fj,)v(ij)dfji (1.22) 



where the resolvent kernel 



is determined by 



B(A ,, ):= EM. (1 . 23) 

A — [i 



F T = (I-K)" 1 ^ = (I + R)f T (1.24a) 
G = g(I-K)- 1 =g(I + R), (1.24b) 



•5 



J. Harnad and A. Its 



where (1 — K) 1 in (1.24a) is understood as acting to the right, while in (1.24b) it acts to the 
left. These quantities similarly satisfy the nonsingularity condition 

F T (A)G(A)=0, (1.26) 

with the diagonal value of the resolvent kernel given by 

R(X, A) = F T '(A)G(A) = -F T (A)G'(A). (1.27) 

The basic deformation formula relating the Fredholm determinant to the resolvent operator 
is 

dlndet(I - K) = -tr ((1 - K)" 1 ^) = -tr ((1 + R)dK) , (1.28) 

where d signifies the differential with respect to any auxiliary parameters on which K may 
depend. 

The determination of R is equivalent to the solution of an associated matrix Riemann- 
Hilbert problem (cf. [IIKS]). Define the r x r invertible matrix H(X) as the following rank-p 
perturbation of the identity matrix I r 

H(X) =I r + 2nif (A)g T (A) = exp 2mf (A)g T (A) (1.29) 

(where the exponential form holds here because the nonsingularity condition (1.19) implies 
that f (A)g T (A) is nilpotent.) The relevant matrix Riemann-Hilbert problem consists of finding 
a nonsingular r x r matrix valued function x(A) that is analytic on the complement of T, 
extends to A = oo off T, with asymptotic form 

X (A)~I r + 0(A- 1 ) (1.30) 

for A — > oo, and has cut discontinuities across V given by 

X-(A)= x+ (A)tf(A), Aer, (1.31) 

where x+(A) and X-(A) are the limiting values of x(A) as V is approached from the left and 
the right, respectively. From the Cauchy integral representation we then have 

X (A)=I r -/ * WfW f ( 'V (1.32a) 
Jr — a 

x-'(A) = i r+ [ 'MfMf (1 . 32b) 
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where the choice of branch x±(f J >) m the integrand is immaterial, in view of (1.19). It follows 
that the p x r matrix-valued functions F(A), G(A) defining the resolvent kernel are obtained 
by applying the matrices x(A), (x T ) _1 (A) to f(A) and g(A), respectively 

F(A) = X (A)f (A), (1.33a) 
G(A) = (x T )" 1 (A)g(A). (1.33b) 

The matrices %(A) , x _1 (A) are thus determined on the complement of V by the integral 
formulae 

*(A)=Ir+ f ^^M dn, (1.34a) 
Jr a — /x 

x - 1(A)=Ir _ /^^IM^, (1 . 3 4b) 
Jr a — n 



with appropriate limiting values as A — >• V. It follows that if we express the asymptotic form 
of x(A) and x _1 (A) for large A as 

oo 

X (A)~I r + ]Tg, (1.35a) 

oo 
J'=l 

the aymptotic coefficients Xj > Xj are given by 

Xj =^A J - 1 F(A)g T (A)rfA, (1.36a) 
Xj = -J X j ~ 1 f(X)G T (X)d\. (1.36b) 

These must of course also satisfy the algebraic recursion relations 

i-i 

Xj = ~Xj -^Xj-kXk , 3 > 1, Xi = -Xi- (1-37) 
fc=i 

In the next section, we specialize to certain particular choices of the functions f(A), g(A) 
determining the Fredholm kernel, and support curves T, for which the computation of the 
Fredholm determinant can be reduced to the solution of an associated isomonodromic defor- 
mation problem. The resulting systems, together with the generalization discussed in Section 
4, extend those obtained in [IIKS], [P] via the Riemann-Hilbert method and in [JMMS], 
[HTW], [TW1-TW3] by other related methods. 
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2. Special Integrable Kernels and Isomonodromic Deformations 

In the following, we choose n = 2m to be even, and let {aj}j=i... n be a set of distinct 
points along T, ordered according to its orientation, and let {r j } J= i v .. m be the ordered se- 
quence of connected segments of V such that Fj has endpoints (a^j-i, «2j)- Let 

B = diag(/3 1 ,- (2.1) 

be a diagonal r x r matrix with distinct eigenvalues {{3 a }a=i---r and {fj, gj}j=i - m a set of 
fixed r x p rectangular matrices satisfying 

gjf fc = 0, Vj,k. (2.2) 

Define 

m 

f (A):=£f^(A), (2.3a) 
j'=i 

m 

g (A):=£g^(A), (2.3b) 
i=i 

where {6 , J (A)} J= i... m denote the characteristic functions of the curve segments {Tj}, and 

f(A) := * (A)f (A), (2.4a) 

g(A):=(^(A))- 1 g (A), (2.4b) 

where 

*o(A) := e XB (2.5) 

is viewed as a "vacuum" wave function. 

Let x±(A) be the solution to the matrix Riemann-Hilbert problem (1.30), (1.31) corre- 
sponding to the choice (2.4a), (2.4b) for the functions f(A),g(A). Define the invertible r x r 
matrix valued function 

*(A):= X (A)*o(A), (2.6) 
with limiting values \I/± on either side of the segments of F 

*±(A) :=X±(A)*o(A). (2.7) 

Then ^± satisfy the discontinuity conditions 

*_(A) = *+(A)# (A) (2.8) 



8 



Integrable Fredholm Operators and Dual Isomonodromic Deformations 



across these segments, where 

H (X) := I r + 2nih (X) = exp2nih (X), (2.9) 

and 

n 

h (X):=J2^sje j (X) (2.10) 

3 = 1 

has rank < p : is piecewise contant along the segments Tj and, because of (2.2), satisfies 

ho(X)ho(fi) = 0, VA,/x. (2.11) 

In terms of these quantities, H(X) is given by 

H(X) = y (X)H (X)^\X). (2.12) 

From (1.33a), (1.33b) and (2.6), it follows that the matrix-valued functions F(A), G(A) de- 
termining the resolvent kernel are given by 

F(A) = *(A)f (A), (2.13a) 
G(A) = (vl> T (A))" 1 go(A), (2.13b) 

where, again, the choice of *±(A) on the RHS is immaterial because of the orthogonality 
conditions (2.2). 

2a. Isomonodromic families of operators: deformation of the endpoints. 

In this section, we construct an isomonodromic family of operators of the form (1.1a), 
(1.1b) having the "dressed" wave function \I/(A) defined in (2.6) as kernel. This follows by 
examining the analytic properties of \I> (A) in a neighborhood of the singular points {A = 
ttj}j=i...n and asymptotically as A — > oo. Specifically, define {-/Vj} J= i... n by 

Nj := -FjGj, (2.14) 

where 

F j := lim F(A) G 3 := lim G(A), (2.15) 

A — >ctj X — >aj 

with the limit A — > «j taken inside the segment r^. We then have the following result. 
Theorem 2.1: The wave function \I> ( A) defined by (2.6) satisfies the equations 

^-(^t^-)*^ (2.16a) 



dX \ ^ X — a 



dctj X — ctj 



* = 0, (2.16b) 



9 



J. Harnad and A. Its 



with the Nj 's given by (2.14). This implies the commutativity 

[D x ,V aj }=0, [V ai ,V a] ]=0, i,j = l,...n (2.17) 

of the operators 



OX A — a* dX 

j=l 3 



V aj :=^- + -^, (2.18b) 
ocxj X — aj 



and hence the invariance of the monodromy data of the operator V\ under changes in the 
parameters {<x,} . 

Proof: The discontinuity conditions (2.8) may be written 

m 

¥+(A) " = -27rz* + (A)^f jg J^(A). (2.19) 

From analyticity, and the fact that the fjgj satisfy 

fjgjffcgfc = 0, Vj,/c, (2.20) 

which follows from (2.2), we conclude that, in a neighborhood of Tj, *(A) can be written 

f s T 



A - a 2 - x 
A - a 2 j-i 



tt(A) = V>,-(A)( /' ) , (2.21) 



where "!/>■/ (A) is locally holomorphic and invertible. It follows that the logarithmic derivative 
^a^ 1 is meromorphic for finite A with simple poles at the points {A = aj} J= i... n , and the 
residue matrices are given by 

:= (-ly^^ia^^^-^a,). (2.22) 

The local form (2.21) impies that these Nj's are given by (2.14), with Fj, Gj defined by (2.15). 
Note that, by eq. (1.26), we also have 

FjGj = 0. (2.23) 

From the asymptotic form (1.30) and the equation (2.6), we see that, for large |A|, ^/a^ -1 
is of the form 

^a^" 1 ~ 5 + C(A _1 ), (2.24) 
10 
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and hence, by Liouville's theorem, \1> satisfies (2.16a). 

From the local form (2.21) of ^ and the fact that \1/ is analytic away from T, we see that 
is holomorphic away from aij and that in a neighborhood of ctj, the sum 



q,- 1 + 



dctj A — OLj 



(2.25) 



is analytic. From the asymptotic form (1.30) and (2.6) it also follows that J^-\I/ _1 vanishes at 
A = oo and hence, again by Liouville's theorem, the expression in (2.25) vanishes identically, 
hence proving that \1/ satisfies (2.16b). 

The fact that the invertible matrix \1/ satisfies both (2.16a) and (2.16b) implies the com- 
mutativity conditions (2.17) for the operators U\, V a . defined in (2.18a) and (2.18b). These 
deformation equations for the operator T>\ imply the invariance of its monodromy data both 
at the finite regular singular points {&j}j=i... n an d at the irregular sigular point A = oo (cf. 
[JMMS, JMU]). Q.E.D. 

Equivalently, we may consider the equations for the matrices Fj, Gj implied by taking 
the limits A — > a k in (2.16a), (2.16b). 

Corollary 2.2: The matrices Fj, Gj satisfy the differential system 



da k 

d_G^ 

da k 

dG^ 

dati 



V 



N k 



\ 



fc= i a 3 - a k 



J ■ 



OLj - OL k 



( 



n N T 



\ 



Gj. 



(2.26a) 
(2.26b) 
(2.26c) 
(2.26d) 



Proof: This follows from the definitions (2.13a), (2.13b), (2.15) and equations (2.16a), 
(2.16b) of the theorem. 

Remark. Let N = np and define the pair of rectangular N x r matrices F, G from the p x r 
blocks {Fj, Gj} 

/Ff\ /Gfx 



F := 



G 



Gl 



\G T J 



(2.27) 
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Then 

AT. 

Af(X) =B + J2— l — = B + F T (A- AI)" 1 ^ (2.28) 

3=1 aj 

as in (1.10), where A is the diagonal N x N matrix with eigenvalues {otj}j=i... n all of multi- 
plicity p, and the equations (2.26a)- (2.26d) are just the Hamiltonian equations (1.11a). 

2b. Deformations of exponents. 

We now turn to the dependence of the wave function \I> (A) on the parameters {P a }a=i,...r- 
By examining the analytic structure of the logarithmic derivatives J^- near T and asymptot- 
ically, we can similarly derive differential equations for \1/ with respect to these parameters, 
and hence deduce the independence of the monodromy of T>\ of the parameter values. 

Theorem 2.3: The wave function *(A) satisfies the equations 

T> Pa V = 0, a = l,...r, (2.29) 
where the operators {^ , / 3 a } a =i,...r are defined by 

8 *E a N,) E b + E b AT,-) E a 

with the Ni's given by (2.14), (2.15) and E a given by (1.4). This implies the commutativity 
conditions 

P>A, Z>/?J = 0, [X>^,X>^]=0, o,6=l,...r, (2.31) 
and hence invariance of the monodromy data of V\ under the deformations parameterized by 

\Pa}a=l,...r- 

Proof: Since the matrix H (X) entering the discontinuity conditions (2.8) is independent 
of the parameters (3 a , we see that the logarithmic derivatives of ^± with respect to these 
parameters have no discontinuity across V 

^ = Jt r ~ - (2 ' 32) 

From the local structure (2.21), it follows that J^- 1 !/ -1 is in fact holomorphic in a neighborhood 
of T, and hence throughout the complex A plane. Furthermore, from the asymptotic form 
(1.35a) and (2.6), it follows that the difference 

^*-i-(AK+[xi, E a }) (2.33) 
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vanishes in the limit A — > oo. Therefore, again by Liouville's theorem this difference, being 
globally holomorphic, vanishes identically 

— -(XEa + bb, E a ])V = 0. (2.34) 

Since \1/ is invertible, the compatibility of (2.34) and (2.16a) implies 



|- fi -EA> ^--^ a -[xi, E a ] 

o\ A — a, c/,y a 



= 0, a=l,...r. (2.35) 



From the asymptotic form of this equation, it follows that the off diagonal terms of xi are 
given by 

n 

(&-&)(xi)«6 = £(^)«6> ( 2 -36) 

and hence 

[xuE a ] bc = (5 ab - 5 ac ) K Jbc (2.37) 

Pb — Pc 

Equation (2.34) is therefore just the condition that $ be in the joint kernel of the operators 
defined in (2.30) and (2.35) are the commutativity conditions implying that the monodromy 
of V\ is preserved under the deformations parameterized by {(3 a } a =i,...r- Q.E.D. 
Remark. The compatibility of eq. (2.34) with (2.16b) also implies the commutativity condi- 
tions 

[V ai , V Pa }=0, i = l,...n, a=l,...r (2.38) 

Finally, we may again equivalently consider the system of equations satisfied by the 
quantities {Fj, Gj}j=i~-n, which follow from (2.16a) and (2.29). 
Corollary 2.4: The matrices Fj, Gj satisfy the differential system 

dF ■ 

= {a 3 E a + [xi, E a ]) F jt (2.39a) 

0Pa 

8C 

1 £ = (-a j E a +[xT, E a ])Gj, (2.39b) 

with the term [xi, E a ] given by (2.37), which are exactly the Hamiltonian equations (1.11b) 
under the identihcations (2.27). 

2c. The Fredholm determinant. 

We now proceed in a standard way to compute the logarithmic differential (1.28) with 
respect to the parameters {aj, /? a }i<j<n, i<a<r- The result is summarized by the formula 
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Theorem 2.5: 

n r 

dlndet(I-K) = u = H kda k + K adf3 a , (2.40) 

k=l a=l 

where the individual factors may be expressed 

gt= 8 lndet(I-K) = ^trm, (2 . 41a) 

Oak ~~f OL-OLj 

Ka _ gj^K) _ t aj W J + ± iEU y a. ^ . ( , 41b) 

"Pa j = 1 6=1 Pa Pb 



Proof: To derive (2.41a) we use the equation 



trTC n 

lndet(I-K) = trln(I-K) = . (2.42) 

n=l 71 



From the formula 



trK n = tr 



rc*2j 1 r a 2j n 

J2 / •••/ K(\ 1 ,\ 2 )K(\ 2 ,\ 3 )...K(\ n ,\ 1 )d\ 1 d\ 2 ...d\ n , (2.43) 



we conclude that 

^ = (-l) k ntvK n (a k , « fc ), (2.44) 

where -fC n (A, /x) stands for the kernel of K n and "tr" in the r.h.s. is the matrix trace. Equations 
(2.42) and (2.44) imply (cf. [TW]) that 

91nde ' (1 ~ K) = (-l) fc+1 tr £ JT» (a fc , a fc ) = (-l) fe+1 tri?(a fc , « fc ), (2.45) 

where the second equality follows from the Neumann expansion of the resolvent operator R. 
From equations (1.27), (2.13b), and (2.16a) it follows that 

/ n N T \ 
R(X, A) = F T (A) f B + J2 XZ^- ) °( A )- ( 2 - 46 ) 

Passing in the last equation to the limit A — > a k and taking into account (2.15) and (2.23), 
we conclude that 



R(a k ,a k ) = (-l) k Fl 



V 



3 r[ a k - aj 



G fc , (2.47) 
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which, because of (2.14), implies that 

( 

tiR(a kl a k ) = (-l) k+1 

\ 



and hence (2.41a). 

To prove (2.41b) we use (1.28), which implies 



From (2.4a) it follows that 



tf = - tr ( (1 -K)-'!f). (2-49) 

= f T (A)£ agM . (2.50) 
Equations (2.49) and (2.50) imply that 

K a = -tv J^F T (/j)E a g(n)d^ (2.51) 

which in view of (1.36a) means that 

K a = -tv(E aX i) = -(Xl)aa- (2.52) 
From the asymptotic form of the equation (2.16a), it follows (cf. (2.36)) that 

n 

[ Xl ,B} = J2Nj, (2.53a) 
j'=i 

n n 

Xi = [X2,B] -J2 N m -Y. a i N i- ( 2 - 53b ) 

3=1 3=1 

The last two formulae determine the diagonal elements of xi as 

(Xl)aa = — a j(Nj)aa ~ " 7T— 5 , (2-54) 



j=l 5=1 



k)bo 



from which follows the expression (2.41b) for K a . Q.E.D. 

2d. Monodromy Data and the Fredholm Determinant as a r— Function. 

Let 

15 
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be an arbitrary rxr linear system having the same singularity structure as the system (2.16a). 
According to [JMU] , the r-function associated to a given solution 

Aj = Aj(ai, ...,a n ;/?i, ...,(3 r ) (2.56) 

of the corresponding isomonodromic deformation equations is defined via the differential form 
(see formula (5.17) in [JMU]): 

1 da — da r n 

lo =-VtrA j A fc -^ -+EE^'W a ^«) 

2 «,■ — at, ^-^ ^-^ 

j^k 3 K o=lj=l 

2.57 

2 a^b j=l j=l Pa Pb 

More precisely, it is shown in [JMU] (and also follows from the Hamiltonian structure of 
the deformation equations [H]) that the form (2.57) is closed, allowing one to introduce the 
r-function, defined up to multiplication by local contants by the equation: 

dlnr = u. (2.58) 

The monodromy data corresponding to system (2.55) consist of the set of n+1 monodromy 
matrices, 

Mj, j = l,...,n, Moo, (2.59) 

M 00 M n M n _ 1 ...Mi = I r , (2.60) 

associated to the singular points, 

aj, j = l,...,n. = oo, (2.61) 
augmented by the two Stokes matrices, 

Si, S 2 (2.62) 

related to 0^, which is the only irregular singular point of (2.55). The Riemann-Hilbert 
problem (2.8) provides a special solution of the isomonodromy equations, 

A j (a 1 ,...,a n ;P 1 ,...,p r ) := Nj, (2.63) 

where Nj is defined in (2.14). This solution is characterized by the following choice of the 
monodromy data (cf. (2.21)): 

M 2j = M^.j = I, 27rif,-gJ = exp(-2 7 rzf J gJ), j = 1, m, n = 2m, (2.64a) 
= Si = S 2 = Ir- (2.64b) 
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Comparing now the equations (2.40) and (2.57) (cf. also eqs. (1.6)— (1.8)), we conclude 
that the r-function evaluated on the solution (2.63) of the isomonodromy equations is given 
by the Fredholm determinant, 

r(ai,...,a n ;/? 1 ...&) = det(I-K) (2.65) 

of the integrable Fredholm operator K characterized by the choice (2.4a), (2.4b) for the func- 
tions f(A),g(A) in (1.18). It is also worth noticing, that the equation (2.65) follows directly 
from the equations (1.6)-(1.8). 

We now show that this Fredholm determinant is also a r-function in the sense of Segal- 
Wilson [SW] and Sato [Sa]. For this, we assume V is a simple, closed curve containing the 
origin of the A-plane in its interior. (In [SW], V is taken as the unit circle, but this is of 
no significance to the general formulation.) Following [SW], for any positive integer k, let 
Ji k = L 2 (T, C fc ) be the Hilbert space of square integrable C fc -valued functions on T, and let 

n k = n k + + n k _ (2.66) 

be its decomposition as the direct sum of subspaces H+, Tt 1 ! consisting, respectively, of el- 
ements admitting a holomorphic extension to the interior (+) or the exterior (— ) of T, the 
latter normalized to vanish as A — > oo. Now, take k = r, and let W C 7i r be a subspace that 
is the graph of a compact linear operator 

h w :H r + ^H r _. (2.67) 

Let 

7 (t) :H r -^H r , 7 G 1^ (2.68) 

denote the action of an abelian subgroup F r C <5l(r) + of the loop group <3l(r)+ of nonsingular 
r x r matrix valued functions on V admitting a holomorphic nonsingular extension to the 
interior of T, acting by left multiplication on H r . Let W(t) be the image of W under the 
action of 7(t), and choose (in the sense of [SW]) an admissible basis, in which W is represented 
by the (2oo) x oo matrix 




(2.69) 



where w + , u>_ represent maps 

w+ : H r + -> H r + , W-:H r + ^ H r _ (2.70) 
the summed images of which give W. In this notation, we have 

W- = hw ° w+ (2-71) 
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The action of 7(t) can then be expressed in matricial form as 
where the maps 



w. 



(2.72) 



a(t) 
bit) 
d(t) 



H r + 



n 



+ 



(2.73a) 
(2.73b) 
(2.73c) 



are defined by restriction of 7(t) to H± composed with the orthogonal projection maps 



P±:H r -> W± 



(2.74) 



a '■= P+ ° "flm, b:=P + o^\ 



d := P- o 7|^r . 



A convenient way to express these projection maps is via the Cauchy operators 



(P ± v)(X) = 



lim 



2ni \±^x J r fj, — X 



2yri 



fi-\ z 



■dfj,, 



(2.75) 



where the expression on the right is an abbreviated notation for the limit, indicating that after 
integration \± approaches A G T from the corresponding side T±. The Segal-Wilson-Sato 
r-function for such a subspace W C H r is then just given by the determinant 



T W {t) := det (1 + 6(t) o h w ° a(t) x ) . 



(2.76) 



For the specific case at hand, we let T be a closed curve passing through the points 
{cKj}i=i... n an d define to be the image of Ti r + under the action of the loop group element 
# _1 (A) where H (X) is defined in (2.9) 



W = W(a u ...a n ) := H-\H r + ) = {v+ - 2irih Q v +1 v+ E H\}. 



(2.77) 



The abelian subgroup T r + C 0l(r) is taken to be the set of values of the "vacuum" wave 
function 

7 (t) := MPU -..&) = e AS , {t a := /?a}a=l,...r- (2.78) 



Then, using (2.75), and the orthogonality condition (2.11), we conclude that 

ho(fj,)v+(n) 



18 



(2.79) 



Integrable Fredholm Operators and Dual Isomonodromic Deformations 



Hence we may express the operator hw as 



h w v + (X) = I hMv ^ d^ (2.80) 
/i — A_ 



r 



while the maps a, b and d are given by 



av+(X) = * (A)v+(A) (2.81a) 

v+eH r + , v_eH r _. 

We then have the following theorem. 
Theorem 2.6: The following two determinants are equal 

det(I - K) = det (I + 6(t) o h w o a(t)" 1 ) , (2.82) 

and hence the r-function t(cki, . . . a n ; , . . . (5 r ) evaluated in (2.65) is the same as the Sato- 
Segal-Wilson r-function Tw(t) defined in (2.76). 

Proof: Composing the inverse a -1 of the map defined in (2.81a) with hw and b as defined 
in (2.80), (2.81b) gives the map b o hw ° « _1 : ri\ — > ^+ as 

» „„, fl - 1(o+)(A) = _1 jf jf MMM^ ( , 83) 

Inverting the order of integration, evaluating the residues at v = fj, and v = X+ and using 
(2.10) gives 

bohwo a-\v + )(X) = -l Wg r (/*K(/*) d/i + / ^(W^KM ^ (284) 

Jr — X + J r [i — A + 

with f(/i), g(/i) defined by (2.4a), (2.4b). Now note that K : W — > W can be written in the 
form 

K = Af o B g , (2.85) 

where B g : W — > 7i!j_ is defined by 

B g 0(A) = / g( f )0M ^, G (2.86) 
Jr X — fx 
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while Af : H r + — > W is defined by 

A f t»+(A)=f T (A) V+ (A), v+eH r + . (2.87) 
On the other hand, from (2.84), we see that b o hw ° « _1 can be expressed as 

bohwoa' 1 =Bf oA g (2.88) 

where 

Bf = B f - a o B fo (2.89) 

with B f and B fo defined as in (2.86), with g replaced by f and f , repectively (cf. (2.4a), 
(2.3a)), while A g is defined as in (2.87), with f replaced by g. 

The Weinstein-Aronszajn identity implies that 

det(I + Bf o A g ) = det(I + A g oB f ) = det(I + A g o B f ), (2.90) 

where the second equality follows from that fact that, due to the orthogonality conditions 
(2.2), we have 

A g o a o B fo = A g0 o B fo = 0. (2.91) 

Therefore 

Tw (t) = det(I + b(t) o h w o a(t)" 1 ) = det(I - K T ), (2.92) 

where 

K T :=-A g oBf. (2.93) 

Comparing (2.85) with (2.93), we see that K T is really just the transpose of K, and hence the 
determinants in (2.82) and (2.92) are the same. Q.E.D. 
Remark. An identity related to (2.90); namely, the equation 

(I + Bf o Ag)" 1 + Bf (I + Ag o B f )" 1 Ag = I, 

lies, in fact, behind the basic relationship (1.33a), (1.33b) between the Riemann-Hilbert prob- 
lem and integrable Fredholm kernels (see [DIZ], proof of Lemma 2.12). 

3. Duality 

3a. The Dual Isomonodromic System. 

In the following, choose r to be even, r = 2s, and let V be a piecewise continuous oriented 
curve in the complex z-plane consisting of the ordered sequence {r a } a= i j ... s of connected com- 
ponents whose endpoints are given by the consecutive pairs {faa-i, P2a}a=i....s of eigenvalues 
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of the qr x qr matrix B = diag(/3i, . . ./32s), where each eigenvalue j3 a has multiplicity q, with 
q < n. Choose a set {f a , g a } a =i,...s of fixed rectangular nx q matrices satisfying 

fjg b = 0, a,b=l,...,s (3.1) 

Let 

s 

f(z) = e zA J2Wz), (3.2a) 

a=l 

s 

z{z) = e- zA Y J ZaO{z), (3.2b) 

a=l 

where 9 a is now the characteristic function, along T, of the curve segment T a and A is the 
diagonal n x n matrix 

A = diag(o:i, • • -a n ). (3.3) 

Let x{ z ) De the nonsingular n x n matrix valued function satisfying the conditions that it be 
analytic in the compactified complex z-plane in the complement of T, with asymptotic form 

X(z) = l n + 0(z- 1 ), (3.4) 

having cut discontinuities across V given by 

X-(z) = x+(z)H(z), zef, (3.5) 

where x+( z ), X-( z ) are the limits as V is approached from the left and right, respectively, and 
H(z) is the invertible n x n matrix function defined along V by 

H(z) = I n + 2nii(z)g T (z) = exp27rif(z)g T (z). (3.6) 

As above, we define the q x q matrix Fredholm integral operator K acting on C^-valued 
functions v on T by 

K(v)(;z) = J^K(z,w)v(w)dw, (3.7) 



where 



K(z,w) = iT{zMw) . (3.8) 
z — w 



The resolvent operator 

R=(I-K)" 1 K (3.9) 
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may again be expressed as 

R(v)(;z) = J^R{z,w)w{w), (3.10) 

where the kernel _ _ 

~R M =* T{Z)&(W) (3.11) 
z — w 



is determined by 



F T = (I — K)" 1 ^ = (I + R)f T = f T x T , (3.12a) 
G = g(I - K)" 1 = g(I + R) = (x T )" 1 g. (3.12b) 



As above, x{ z )iX 1 (. z ) have the asymptotic expansions 



xW = i» + Ef' (3 - 13a) 

oo ~ 



where 



satisfy 



Xi = J^z^ 1 F(z)g T (z)dz, (3.14a) 

Xj = - J^z j - 1 f(z)G T (X)dz. (3.14b) 

Xj = -Xj - 5^ Xi-fcXfc > J > 1) Xi = -Xi- (3-15) 
fe=i 

Again, we may define the n x q rectangular matrices 

F a := lim F(z) G a := (-l) a lim G(z), a = l,...r, (3.16) 

where the limits are taken inside the curve segments {F a }a=i,.. s , the n x n matrices 

M a :=-¥ a G T al a = l,...r. (3.17) 

and the pair of n x rq matrices F, G formed from the blocks {F a , G a } 

F:=(F! F 2 • • F r ), (3.18a) 

G:=(Gi G 2 • • G r ). (3.18b) 
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Then 



M(z):=A + J2 

o=l Z Pa 



-1P1T 



A + F(B — zl)~ G 



as in (1.12a), (1.12b). 

Now define, similarly to the previous section, 

:=x{z)e zA , 

with limiting values on either side of the segments of F 



*±(*) -=X±(z)e zA . 



(3.19) 



(3.20) 



(3.21) 



It then follows from the same arguments as in the preceding section that \j> simultaneously 
satisfies 



V z ^ = 0, 
V 0a $ = O, a=l,...r, 
P Qi * = 0, i = l,...n, 



(3.22a) 
(3.22b) 
(3.22c) 



where 



V, 



V 



d_ 

dz 

d 



-M(z) 



d_ 

dz 



- A -E 



Pa 



+ 



M n 



d(3 a z-/3 a 



d 



M a 



da, 



-zEi-Y, 



3=1 

3*i 



Ei (Ea = i M «) Ei ± E i S=i M ") E ' 



(3.23a) 
(3.23b) 
(3.23c) 



with {-Ei}i=i,...n the elementary nxn matrices with entry 1 in the (ii) position. Consequently, 
all these operators commute and the monodromy of the parametric family of operators V z 
is invariant under the changes in the deformation parameters {aii, P a }i=i...n, a—i,...r- 

The 

corresponding dynamical equations, expressed in terms of the matrices {F a , G a }a=i,...r are 



dF a 

dF a 

d(3 a 



F a , a ^ b, 



( 



\ 



(3.24a) 
(3.24b) 
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<9G 



dFa 

don 

dGa 

dai 



V 



6^o 



/ 



1 Pa- Pb 

a 

(0 a Ei + \xi, Ei])F a , 
(-P a E i+ [xf, Ei])G a , 



G a , 



(3.24c) 
(3.24d) 

(3.24e) 
(3.24f) 



with the term [xi, Ej\ given by 



[Xi,Ei] jk = (Sij - 5 lk ) 



(Ea=l M a ) 



jk 



a.j - a k 



(3.25) 



These are again just the Hamiltonian equations (1.11a), (1.11b) under the identifications 
(3.18a), (3.18b), provided the pair (F, G) is identified with (F, G). 

Finally, we may compute the logarithmic differential of the Fredholm determinant of K 
to obtain, as in Section 2c, 



rflndet(l - K) = u = ^ H j da j + ^2 K °- d Pa, 

3=1 



(3.26) 



a=l 



where 



VP a b=1 Pa Pb 



(3.27a) 



gjn det(I - K) = ^ I (SU (ELi ^ 

act,- a, — «fc 

J a=l fc=l J K 



3b. Duality Theorem. 

In order to relate the above isomonodromic deformation system to the one in Section 2 we 
henceforth make the following restrictions. First of all, we consider only the case p = q = 1, so 
K, K are both scalar Fredholm integral operators. The rectangular matrices f (A), g(A), fj, gj 
all become r-component column vectors while f (z), g(z), f a , g a are n-component column vec- 
tors. Furthermore, the fj's and f 's are chosen to be of the special form 

f, = e, f a = e, j = l,...m, a=l,...s, (3.28) 
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where e e C r , e e C n are the column vectors with all components equal to 1, while the 
components of the column vectors gj and g a are chosen such that 

(gj)2a = -(g J ) 2 a-l = (ga)2j = -(ga^j-l == C ja (3.29) 

for some constant m x s matrix with elements {cj a }j=i,...m,a=i...s- 

We then have the following result, which shows that the isomonodromic deformation 
system (3.23a)- (3.23c) is, in fact, the dual system to the one defined in Section 1 by eqs. 
(2.18a), (2.18b), (2.30). 

Theorem 3.1: K is the Fourier-Laplace transform of K along the curves V and V and 

F = F, G = G, (3.30a) 
u = u (3.30b) 

Proof: On the product T x T, define the following locally constant function 

m s 

K{\z) :=^^c^(A)0 a (*). (3.31) 

3 = 1 a=l 

Taking the Fourier-Laplace transform with respect to the variables z and A along the curves 
T and T, respectively, gives the two Fredholm kernels 

K(X,fi)= [K(u,z)e zix -^dz= fT [ X ^\ (3.32a) 
Jr A — [i 

K(w,z)= [ K(u,z)e^ w - z) du.= fT Hg(^ (3.32b) 
Jr w — z 

where f, g, f, g are defined by (2.4a), (2.4b) (3.2a), (3.2b) with gj f a , g a given by (3.28), 
(3.29). 

To prove (3.30a), we use the Neumann expansion of the resolvent R in (1.24a) to get the 
following convergent series for the components of F 

F a ( aj ) =Ua 3 ) + J K(a j ,X 1 )f a (X 1 )dX 1 + J^J^ K(a 3 , X 2 )K(X 2 , X 1 %(X 1 )dX 1 dX 2 + ■■■ 

J J„K(X U z^e^-^e^dzxdXx 

III lk(X2,Z2)K(X 1 ,z 1 )e Z2{ ^- X2) e Zl{X2 - Xl) e^ Xl dz 1 dz 2 dX 1 dX2 + --- 
J r «/ r r r 

=e a ^ a + jf, J k(X 1 ,z 1 )e Xl( - (3a - Zl) e Zia >dX 1 dz 1 

+ J~J~J J k(X 2 ,z 2 )k(X 1 ,z 1 )e X2{zi - Z2) e XliPa - Zl) e a i Z2 dX 1 dX 2 dz 1 dz 2 + --- 

=%(Pa) + J^K((3 a ,z 1 )%(z 1 )dz 1 + J^J^k(f3 a ,z 1 )k(z 1 ,z 2 )%(z 2 )dz 1 dz 2 + --- 
=F 3 (p a ), (3.33) 
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where the Neumann expansion of R in (3.12a) is used in the last line. Similarly, using the 
Neumann expansions of R and R in (1.24b), (3.12b) and taking into account that 



X> fc[a ji]0 fc (A) )K( aj ,z) = k(X,f3 a ) [J2 c [i±±]M*)) (3-34) 

fc=l / \6=1 / 

a = l,...r, j = l,...n, 

gives the equality 

G a {a j ) = G j {p a ). (3.35) 

The equality (3.30b) follows directly from the fact that K and K are related by the Fourier- 
Laplace transforms (3.32a), (3.32b). It may also be seen from the equalities (3.30a) in view 
of (2.28), (3.19) and the general duality relations mentioned in the introduction and proved 
in [H]. Q.E.D. 

4. Reductions and Generalizations 

4a. Symplectic reduction. 

In applications, the structure of the Fredholm integral operator K and the associated 
Riemann-Hilbert problem data may be such that the generic systems appearing in sections 
2a, 2b are reduced to systems having fewer independent variables. An important example of 
this involves reduction from the group (5 [(2s) to the symplectic subgroup &p(2s). Let J be 
the symplectic 2s x 2s matrix having block form 



. l s o ; ' 

and suppose that the matrix B of section 2 satisfies the relations 

B T J + JB = 0, (4.2) 

implying it is in the symplectic subalgebra 5p(2s) C $j[(2s). Since B is diagonal, with eigen- 
values {(3 a } a =i,...2s, this just means that 

(3 a+s = -(3 a , a = l,...s, (4.3) 

so B has the form 

n = ( 

-By 

where Bi = diag(/3i, . . . f3 s ). This implies that the vacuum wave function ^o(A) satisfies 

*o (A) J* (A) = J, (4.5) 



B=(i 1 ° ) (4.4) 
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and hence takes values in <5p(2s). We also require that the r x p matrices {fj,gj} entering in 
(2.4a), (2.4b) satisfy the relations 

g J f/j+ Jf,-gj = 0, (4.6) 

implying that the group element H Q (X) defined in (2.9), (2.10) satisfies 

HZ(\)JHo(\) = J, (4.7) 

and hence also is in &p(2s). It then follows that the dressed wave function \I/(A) obtained 
by solution of the associated Riemann-Hilbert problem also takes values in &p(2s), and the 
residue matrices Nj entering in the definition of the operator T>\ have values in sp(2s). 

More explicitly, assuming the (fj, gj)'s all have rank p, eq. (4.6) implies that there exist 
symmetric, invertible matrices Sj G <3l(p) such that 

fj = JgjSj, sj = Sj, j = l,... s. (4.8) 

The orthogonality conditions (2.2) then reduce to the symplectic isotropy conditions 

gjjg fc = 0, j = l,...s. (4.9) 

Since the Sj's are invertible and symmetric, we may take their square roots (determined 
up to a Z£ ambiguity) and, without changing the content of the Riemann-Hilbert problem, 
redefine the fj-'s and g/s by absorbing these square roots through the substitutions 

f;^f^J> Sj^KjsJ*- (4-10) 
This just amounts to setting the s/s equal to the identity element, reducing (4.8) simply to 

fj = J Sj . (4.11) 

Since the solution of the Riemann-Hilbert problem preserves the symplectic reduction, it 
follows from eqs. (2.13a), (2.13b) and (2.15) that the matrices (Fj, Gj) satisfy the relations 

F, = (-lyjGj. (4.12) 

This means that there exist pairs of rectangular p x s matrices (Qj, Pj) such that the Fj's 
and Gj's may be expressed 

Fj = (P„ -Qj) , Gj = (-iye^ (Q„ P,) , j = l...n. (4.13) 
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The residue matrices Nj therefore have the reduced block form 



showing explicitly that they all belong to sp(2s). 

The reduced form of eqs. (2.26a)- (2.26d) then becomes 



^ t= (P J Q?-Q < PDQ t 
«P I= (P i Q?-Q J P?)P t 

= -Q,-Bi - V ±-£<h ^ k) ^\ 4.15c) 

0«j aj-afc 



(4.15d) 



In particular, for p = s = 1, r = 2, these are just the equations of [JMMS], Theorem 
7.5, determining the n-particle correlation functions for an impenetrable bosonic gas (or the 
spectral distribution function for random unitary matrices in the double scaling limit [TW1]). 
Similarly, the reduced form of (2.39a) and (2.39b) is given by 

**i = - ajQjEa + Q 7 Y E « Np Q Eb ± EbNpQE a 

p \ E a N QQ E b + E h N QQ E a (AiR a \ 



6=1 



dPj ^E a N QP E b + E b N QP E a 
— =a jPj E a - P, ^ _ 

>A E a N PP E b + E b N PP E a 

Qj S — — ( ' 

a=l,...s, j = l,...n, 

where 

n n n n 

*QQ' ^Ql'Q,- Npp-.^PjPj, N QP :=J2QjPj, N pq := £ PjQj, (4.17) 

j'=i i =1 j =1 j =1 
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and E a is the elementary s x s matrix with 1 in the aa position. 

Other discrete reductions of the generic Riemann-Hilbert data and the corresponding 
isomonodromic deformation families may similarly be derived along the lines indicated in [H] . 

4b. Higher order poles. 

Using arguments similar to the Zakharov-Shabat dressing method (see e.g. [NZMP]), it is 
possible to generalize the above considerations to a broader class of isomonodromic deforma- 
tion equations associated to matrix Riemann-Hilbert problems of a similar nature, by allowing 
higher order poles at one or more of the singular points. For simplicity, we restrict attention 
to the case where there is just one irregular singular point of arbitrary index located at oo. 
Extension to the case of any number of irregular singular points is quite straightforward. 

Let \I/o(A,t) € <5[(r) be an isomonodromic "vacuum" solution, satisfying 

d^ 



dX 
d^ 



=E/o(A,t)* , (4.18a) 
=V°(A,t)*o, a = l...k (4.18b) 

where {Uo(X, t), V^°(A, t)}j=i...fc are r x r matrix valued polynomials in A whose coefficients 
are functions of the /c-component vector t = (ti,...tk)- The compatibility conditions for 
(4.18a), (4.18b) imply that the generalized monodromy data for the "vacuum" operator 

V° x :=^-U (X,t) (4.19) 

are invariant under the deformations induced by changes of the parameters (t±, . . . £&) satisfying 
these equations. For example, we could consider the simplest case, when the matrices V^°(A) 
are independent of the parameters (ti, . . .tk), and commute amongst themselves 

K(A),H°(A)] =0, a,b=l,...k. (4.20) 

The vacuum wave function then has trivial monodromy and can be normalized as 

*o(A,t) =exp (j>aK (A)j , (4.21) 

so that its values form an abelian group. More general cases, in which the vacuum monodromy 
is not necessarily trivial may also be included (cf. [TW2-3]). 

As in section 2, choose a set {fj,Sj}j=i—m °f fixed r x p rectangular matrices satisfying 
(2.2), and let 0j(X) denote the characteristic function along the curve segment r^. Define 

m 

f(A,t):=*o(A,t)^f^-(A), (4.22a) 
j'=i 

m 

g(A,t) := K(A,t))- 1 ^g J ^(A), (4.22b) 
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as in (2.4a), (2.4b), but with the exponential vacuum wave function \I/o(A) replaced by the 
vacuum solution \I/o(A,t), evaluated along the curve segments Tj (where, in the event that 
there is nontrivial vacuum monodromy, a specific branch is taken). Now define Hq(X) again 
as in (2.9) and H(X,t) as in (1.29) and (2.12), but with the functions f(A), g(A) of eqs. 
(2.4a), (2.4b) replaced by f(A,t), g(A,t) and \I/o(A) replaced by \E f o(A, t). Let x(A) again 
be a solution to the matrix Riemann-Hilbert problem as defined in (1.30), (1.31), with the 
appropriate substitution of H(X,t) for H(X). Then, just as in sections 2a, 2b, by considering 
the local structure near the singularities and the asymptotic behaviour at oo, it is easily shown 
that 

*(A):= X (A)*o(A,t) (4.23) 

satisfies the equations 

N- 

^-C/(A,t)*-^ T -i-* = 0, (4.24a) 



d\ ^— ' A — a-i 



J= i ~i 



0f K(A,t)* = 0, a = l...k, (4.24b) 

N- 

7— + — -^*=0, j = l...n, (4.24c) 
octj X — aj 

where, as before, Nj, Fj, Gj are again given by (2.14), (2.15), while 

U(X,t)= ( X (A)t/ (A,t) x - 1 (A)) + , (4.25a) 
K(A,t) = ( x (A)T/°(A,t) x - 1 (A)) + , a=l,...k, (4.25b) 

where ( ) + denotes projection to the polynomial part in A (cf. [NZMP]). It follows that the 
operators 

r\ Th tit 

V X : = U(X,t)-j:^- (4.26a) 

j = l J 

V ta : =A_y a (A,t), a = l...k, (4.26b) 

3 N- 

T> a j :=^- + ^^, j = l...n (4.26c) 
3 octj X — aj 

all commute, and the generalized monodromy data for T>\ is invariant under the deformations 
induced by changes in the parameters {oij}j = i... n , {t a } a =i...k- 

The r function associated with these deformations may once again be shown to be the 
Fredholm determinant det(I — K) of the integral operator K defined as in (1.17), (1.18), with 
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the appropriate substitutions for f(A), g(A). The coefficients {-Hj}j=i... n {K a } a =i...k of the 
differential form 

n k 

dlndet(I-K) = ^Hidai + ^Kadta (4.27) 

i—l a=l 

may again be viewed as commuting nonautonomous Hamiltonians whose Hamiltonian equa- 
tions determine the deformation conditions corresponding to the commutativity of the oper- 
ators V\, T> ta , T> a .. For the specific case (4.21), where the vacuum wave functions form an 
abelian group, the same interpretation may be given to the corresponding r-function as in 
Theorem 2.6. More generally, if we drop the requirement that (2.72) define an abelian group 
action on the Hilbert space Grassmannian 7i r , and simply view the maps determined by (2.72) 
(2.81a), (2.81b), (2.81c) as defining a set of integral curves for a compatible nonautonomous 
system generated by the solution \I/o(A, t) of the associated vacuum system, the same geomet- 
rical interpretion in terms of determinants of projection operators is still valid. The resulting 
differential operators V\, V z may also be expressed in terms of a Hamiltonian quotienting 
procedure through formulae similar to (1.10), (1.12b), and we may thereby again associate a 
dual system and corresponding dual integral operator K, as was done in section 3. Details of 
this, as well as generalizations involving multiple higher order poles will be given elsewhere. 
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